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In the past two decades, the field theory of disordered 
electronic systems, a matrix version of the nonlinear cr- 
model, has attracted a great deal of theoretical attention 
[Q-^. The cr-model provides the perhaps most rigorous 
and comprehensive theoretical framework by which in- 
teracting and non-interacting disordered electronic sys- 
tems can be explored. Three different, albeit closely re- 
lated versions of this model exist: The supersymmetric 
(SUSY) 0, replica || and dynamic or Keldysh |] a- 
model. Due to different microscopic starting points these 
three formulations cover partly complementary applica- 
tion areas: whereas SUSY is tailor made to the analysis 
of non-perturbative phenonmena, the replica, and in par- 
ticular the Keldysh formalism are strong at problems in- 
volving interactions - about which SUSY has nothing to 
say. This constellation amounts to some theoretical vac- 
uum because sooner or later non-perturbative interaction 
related problems will come into focus. 

An important step towards improving the situation 
was made recently Q (first within the framework of 
SUSY) when it became clear that the large energy asymp- 
totic of Wigner-Dyson (WD) level statistics - the non- 
perturbative prototype problem - may be obtained from 
a careful analysis of the saddle point structure of the 
cr-model (in contrast to a full integration over the field 
manifold which is out of question for the non-SUSY rep- 
resentants). Subsequently these ideas were adopted to 
the replica cr-model |^ , which proved to be useful for the 
solution of simple interacting problems Yet, the ap- 
plication spectrum of the replica theory, excluding non- 
equilibrium phenomena, is significantly narrower than 
that of its Keldysh counterpart. Furthermore, replica 
analyses categorically rely on a formidable analytic con- 
tinuation procedure, which tends to obscure the underly- 
ing phyiscs and essentially complicates practical applica- 
tions. It is therefore desirable to formulate a prescription 
whereby non-perturbative information can be obtained 
from the conceptually more transparent Keldysh formal- 
ism. This is the subject of the present Letter. 

All the theories mentioned above share the normaliza- 
tion condition Z = 1, where Z is the functional parti- 
tion function. It is this universal normalzation property 
which makes the models useful in the analysis of disor- 
der properties The normalization in turn is due to 
a global internal symmetry underlying the model: su- 
persymmetry, replica permutation symmetry or, within 



the Keldysh formulation, a symmetry principle related to 
the causality of the theory. Building on these structures, 
the saddle point analysis of any cr-model must consist 
of three parts: (i) identification of all saddle points and 
their action, (ii) proof that Z = \ holds within the sta- 
tionary phase scheme (iii) saddle point evaluation of any 
specific correlation function. Using WD spectral statis- 
tics as a example we are going to apply this program to 
the Keldysh model of unitary (broken time reversal in- 
variance) symmetry. The generalization to different sym- 
metries and observables will be discussed elsewhere 'm. 

We wish to compute the two-point correla- 
tion function i?(ei,e2) = ((p(ei)p(e2))) = 
AV(27r2)sR((tr{G'(e+)}tr{G(e^)})) of the density of 
states p{e) for a model of free electrons with broken 
time reversal invariance and subject to a random po- 
tential. Here A is the mean level spacing and ((...)) 
stands for a (cumulative) disorder average. Our starting 
point is the low energy effective partition function 
Z = / PQe^^[<3]^ where 

iS[Q] = tr {DdQdQ + AieQ} , (1) 

D is the diffusion constant and L the system size. The 
matrix-valued field Q — {Q'''^/} acts in a product space 
defined through an index I = 1,2 labelling the for- 
ward and backward Keldysh time contour, C' = {< G 
(=Foo,±oo)}, and the energy indices, e, Fourier conju- 
gate to the time variables t on C'. For all what follows 
it will be convenient to discretize the (a priori contin- 
uous) energy variables in units of some small spacing 
5^. Since the application range of the effective action 
is limited to energies |e| < lUiHi this manipula- 

tion implies that the Q's become matrices of finite di- 
mension 2K, where K = 2/{S^t). The internal struc- 
ture of these matrices is defined through (a) the con- 
straint Q^{r) — 1 and (b) Hermiticity, = Q. Finally, 
e — diag(e^,e^) is a 2iir-dimensional diagonal matrix 
where = diag(ei, . . . e„, . . . e^) -i(-l)'0. 

The definition of the low energy effective action (|^) im- 
plies one more important structural element which is not 
explicit in the notation: by definition of the trace 'tr', 
the energy integration/summation over any continuous 
and analytically benign function of the energy indices e 
produces zero (e.g. tr {e} — 0). For lack of better termi- 
nology we will refer to this feature as 'causality'. That 
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the causahty criterion is not exphcit in the definition of 
the trace is not just a matter of notational convenience. 
The crux is that to rigorously retrieve the full causality 
properties of the microscopic Keldysh partition function, 
all energies, e, including energies in excess of the width 
of the spectrum of the microscopic Hamiltonian, have to 
be taken into account. The philosophy behind declaring 
the causality principle to an intrinsic feature of the effec- 
tive action is that S merely represents the low energy 
sector of some larger parent theory, S + S'high- After the 
inclusion of the high energy sector S'high, all energy sum- 
mations could, in principle, be extended to infinity and 
the correct spectral structures of the Keldysh partition 
function would be retrieved. 

This anticipation has been at the root of previous ef- 
fective action formulations of the Keldysh approach, and, 
needless to say, is of relevance for all matters related to 
spectral statistics. We will therefore subdivide our anal- 
ysis into two parts. Taking a pragmatic point of view we 
will first show how the spectral correlation function can 
be obtained from the low energy model, Eq. (|]) , once the 
causality assumption has been made. In a second part 
we will then show how causality can be made manifest, 
on the expense of including large energies. We emphasize 
that part II of the analysis has been included for reasons 
of conceptual completeness; in practical applications, the 
large energy sector will normally not play a role. 

Part I: The simplest way of computing i?(ei, €2) from 
the functional integral over the effective action iS'[Q], 
is to couple the energy vector e to sources. This can 
be done by generalizing e' e'^" = e' + k', where 
the energy diagonal matrices k' = {Snn'^nniKi} and e„j 
are the discrete energies closest to the arguments e;, 
1 — 1,2. It is then a straightforward matter to verify 
that the definition of the Keldysh partition function im- 
plies i?(ei,e2) = Ay{27:^)RedU^^^Z[k]-dis, where 
'dis' stands for the disconnected part of the functional 
average. To keep the presentation simple we will focus in 
the quantitative analysis of this expression on the con- 
tribution of the spatial zero mode Q{r) = Q — const, 
(which ought to reproduce WD statistics). The inclusion 
of the spatially fluctuating modes, which is straightfor- 
ward and does not introduce conceptually new elements, 
will be briefly discussed in the end. 

The key to understanding the structure of the zero 
mode integration lies in the observation that the (source- 
free) effective action iS[Q] — —inA^^ tr {eQ} possesses 
a multitude of 2^^ isolated saddle points. Indeed, 
any of the (energy diagonal |^) configurations A = 
diag (±1,±1,...,±1) solves the stationary phase equa- 
tion of the model, (^(^^[Q] = TrA~^[e,Q] = 0. In what 
follows we are going to show that a Gaussian integration 
around these saddle points produces WD statistics. 

To prepare the integration, consider the contribution 
of any saddle point A with K+p entries -1-1 and K—p en- 
tries — 1 . We first re-order these elements (through some 
global unitary transformation) such that A assumes the 
form A = diag (1, . . . , 1, —1, ... — 1). Next, a set of field 



configurations weakly fiuctuating around A is introduced 
through Q — TAT"^, where the unitary rotation matri- 
ces T = exp ^ _^t^) ; i? is a (K+p) x (if— p)-dimensional 
complex generator matrix, and the block decomposition 
corresponds to the signature of A. 

In principle we should now proceed by expanding the 
action to second order in B and integrate. Fortunately, 
however, there is no need to carry out this program 
for all 2^^ saddle points explicitly. The reason is that 
among the entity of saddle points A, there is one ele- 
ment Aq' = (— 1)'"'"^(5''' that plays a distinguished role. 
What makes Aq special is that, unlike the other A's, its 
structure is compatible with the signature of the imagi- 
nary increments of e. Building on this feature, previous 
analyses of the Keldysh cr-model indeed focused on a per- 
turbative expansion around Aq and did not take the other 
saddle points into account. 

As a warm-up, let us outline how an integration 
around the standard saddle point produces the unit- 
normalization of the source-free partition function Z[k = 
0] = 1. Substituting Q = TAqT"^ into the zero mode 
action and expanding to second order in B we obtain the 
quadratic action 



Integration over B then leads to 
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Zo = const. X e-""^ 'tf{«Ao}^g ■. Fo = JJ (4 - 4^'^^ 

n.n' 

as the contribution of Aq to Z[0]. Due to the 
causality property, tr{eAo} = 0. Similarly, Fq — 
exp{- J2nn' ^'^(''t ~ ^n')} = exp{0} = 1, whcrc the pres- 
ence of the imaginary increments guarantees that the 
branch cut singularity of the logarithm is not touched. 
(We re-emphasize that, at this stage, the causality rule 
has the mere status of a working assumption. In part II 
we will make up for its precise formulation, and show that 
the proper ultraviolet extension of the fiuctuation deter- 
minant fixes the unspecified 'const.' to unity.) Com- 
bining factors, we find Zq = 1. As a corollary we re- 
mark that Zq = 1 implies vanishing of the total contri- 
bution of all non-standard saddle points to Z[0]. Be- 
fore turning to these other saddle points, let us discuss 
the contribution Rq of the standard saddle point to the 
spectral correlation function. A straightforward expan- 
sion of Z[k\ to first order in ki and K2 yields Rq = 

^(Erm' Bn,nBl„^Bl^^,Bn,ri2)B, whcrc stands for 

the Gaussian average over the action (||). Integration 
over B then leads to Rq — l/(2s^), where s — 7rcj+/A 
and u+ = e+ - e," 

We next turn to the discussion of the other saddle 
points. In fact we will focus on just a single non-standard 
saddle point A, namely the one where the signs of the two 
entries corresponding to our reference energies e„j 2 , are 

flipped: A = Aq — 2(5'M„„i + 2S'^Snn2- From the analysis 
of A, the role of all the other saddle points will become 
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clear. The re-ordering needed to bring A into the canoni- 



cal form imphes that the action iS*^ differs from iS)^^ in 
two respects: first, the contribution of the saddle point 
itself ^^[A] = -i7rA~Hr{eA} = i7rA~i2a'+ no longer 
vanishes. Second, in the fluctuation contribution, the 
two energy arguments e„2 and e„j are exchanged. Given 
these structural changes we find it more convenient to 
reverse the order of the evaluation of the functional in- 
tegral: first integrate out fluctuations, then expand in 
the sources. The integration over leads to a fluctu- 
ation factor F similar to Fq above, only that e„j^ and e„2 
are exchanged and coupled to the respective sources ki 2: 
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F = Fn 
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We next add to the products the "missing" factors 
n = TT'1,2 and use that the now unconstrained products 
over n, n' , as well as the factor Fq, give unity (causality). 
This leads to the result 
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where the last equality is valid to leading non-vanishing 
order in the sources, Ki,K2- Notice the proportionality 
of F to KiK2. This implies (a) that the sources in the 
action S[A] may be set to zero (we arc differentiating the 
functional at ki — 0) and (b) that only the non-standard 
saddle point A contributes to the correlation function. 
Indeed, for any other non-standard saddle point one or 
several signatures corresponding to energy arguments 
En 7^ eni.na are changed. These energies are not coupled 
to sources. Repeating the steps outlined above one finds 
that the _B-integral around these saddle points gives zero 
(i.e. the k ^ limit of the factor F above). The same 
argument also shows that the non-standard saddle points 
do not contribute to Z[0]. Differentiating Eq. (^) w.r.t. 
Ki_2 and adding the contribution of the standard saddle 
point, we obtain the well known result 

I — exp(2is) ^ ' ^ ^ 



I 



i?(a>) = --Rc. 
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for the two point correlation function of the zero mode 
theory. We finally mention that the inclusion of spatially 
fluctuating modes into the formalism (a) does not change 
the saddle point structure and (b), after Gaussian inte- 
gration, leads to a renormalization F — > 'D{lo)F , where 

(i?q2)2 
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and q are the quantized momenta associated to fluctua- 
tions in a flnite size system. Combining this with the con- 
tribution of Aq, we reproduce the familiar result |§,||j5| 
for the level statistics of the unitary disordered electron 
gas. Since this result was obtained in a saddle point ap- 
proximation, its validity is restricted to energies uj ^ IS.. 
Indeed, in the opposite limit the fluctuation modes B 
become too light to be treated in the Gaussian approx- 
imation. The fact that the zero mode result Eq. (H) is 



actually exact for any w is a "coincidence" (see however 
Ref. 0), speciflc to the unitary ensemble. 

Fart II: The analysis above crucially relied on the 
causality postulate: summations over continuous func- 
tions of e produce zero. As mentioned above, this fea- 
ture could not be made explicit because the validity of 
the action (|l|) is limited to low energies. In this second 
part of the Letter we are going to construct an energeti- 
cally enlarged formulation whereby large and small e are 
treated on the same footing. This will make the causal- 
ity property manifest. At the same time we will see why 
the pragmatic scheme employed above is sufficient for the 
calculation of low energy observables. 

In order to not unnecessarily complicate the discus- 
sion we will formulate this part of the analysis for an A^- 
dimensional random matrix theory (RMT) Hamiltonian 
H = {-ff^i/} defined through the Gaussian correlation law 
{H^,^) = and {II^,^H^, f^,,) = N-^S^^,S^^,. The advan- 
tage gained is that H has neatly defined universal large 
energy properties, i.e. that we will not need to consider 
the non-universal UV asymptotics of the free-electron 
Hamiltonian underlying Eq. (^. Later on we will argue 
that, as far as the present discussion is concerned, the 
specific modeling of the Hamiltonian is of no relevance. 
Further, to deal with a manifestly [/^-regularized model, 
we compactify our energy variables. This can be done by 
discretizing the temporal Keldysh contours to a lat- 
tice of small spacing 6t. As a result, the energy variables 
e„ G {S„ 2S„ KS,}, where K = 2tt/{6A) > 1. 

The effective Keldysh action for the RMT model 
can be obtained by a straightforward adaptation of 
previous derivations of the a-model for RMT Hamil- 
tonians to the specifics of the Keldysh (T-model. 
As an intermediate result one obtains the partition 
function Z[k] = J exp{iS'[(3, k]} where the ac- 
tion iS[Q, k] = -N (i tr (Q^) _ tr In [z + Q]) with 
z = CT3(5j"i(e-''^35t« _ 1) (cTg acts in the Keldysh 2x2 

/ 1' 

space). At this stage Q = {Q^^,} is a 2ii'-dimensional 
matrix that has been introduced to decouple the H- 
averaged action ||l^ (no constraint = 1 as yet). The 
unusual phase-type appearance of the energy argument 
is due to the time discretization. However, in the limit 
e(5t <C 1, 2 -^■ —ii and we retrieve the standard form of 
the RMT (7- model action [Q. We next subject the action 
to a saddle point analysis (stabilized by the parameter 
N ^ 1). Variation of the action w.r.t. Q yields the 
quadratic equation Q = {z + Q)~^. The 2^^-io\d degen- 
erate energy diagonal set of solutions, A(')(e„) = A^^^ 
is given by 



(7) 



This expression determines the entire spectral structure 
of the model. First, notice that for energies e„ ^ 1 ^ 
5^'^, A(')(e„) = ±1 -I- ©(e^). This means that the solu- 
tions A(')(e„) represent an UV extension of the saddle 
points A discussed in part I. We next ask whether the 
two sign alternatives in Eq. (0) are equivalent or whether 
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the model has a preferred choice. Indeed, the latter is the 
case: for energies e„ » 1 greatly in excess of the width 
of the spectrum, A(')(e„) must approach zero - the free 
Gaussian saddle point of the non-disordered model. This 
condition determines Ao(e,i) — ■^[— z(e„) + <7z{—z^{tn) — 
4)^/^] as the canonical solution. For low energies Ao(e„) 
reduces to the saddle point Aq discussed above. This sad- 
dle point has the important property, X]n[-^o(^")]'^ ~ 
k a positive integer. The outline of the proof is as follows: 
(due to the presence of a finite imaginary increment) the 
summation over e„ is equivalent to an integration of the 
variable w = exp{— Z(5te„} over the complex unit circle. 
It is straightforward to verify that for \w\ > 1, Ao(w) is 
analytic (the branch cut singularity of the square root 
lies inside the unit circle). Further, for \w\ ^ 1, the in- 
tegrand decays as From Cauchy's theorem we 
conclude that the summation gives zero □. Summarizing, 
we have found a complex saddle point structure which ex- 
tends the low energy saddle points A = ±1 discussed in 
part I into the UV regime. The complex structure of the 
theory entails the existence of a 'natural' saddle point 
Aq. Existence and behavior of Aq vitally depend on the 
large energy, e ^ 1, asymptotics of the theory. 

To more explicitly establish contact with the low en- 
ergy regime discussed in part I, we next introduce fluctu- 
ations around the saddle points (^. Dcflning Q = A + P, 
where P is some Hermitian matrix fl^ , and expand- 
ing the action 5 [A + P] to second order in P we ob- 
tain iS'j^^P] = -f tr(p2 + PAPA). This is the UV- 
extension of the low energy action discussed in part I. 
Indeed, substituting the small e asymptotics of A and us- 
ing that for the RMT model, A = tt/A^ (H), we find that 
for energies e -C 1, reduces to the action, Eq. (||). 

One can now step by step repeat the analysis that led 
to the correlation function of part I. The only difference 
is that instead of energy denominators ~ (e^ — e~,) con- 
structions like 1 -I- AQ^^(e„)AQ^^(e„') appear. Due to the 
compact phase- type appearance of the energy arguments 
in Ao(en) all energy summations converge. Further, the 
properties of Aq discussed above imply the vanishing of 
expressions like J2n /(^o(en)) where / may be any ana- 
lytic function. This implements the causality principle. 
In parentheses we note that the detailed execution of this 
program yields a unit normalization of the partition func- 
tion, without any undetermined prefactors. 

The discussion above provokes the obvious question 
whether the conclusions drawn from the high energy 
asymptotics of the action are specific to the random ma- 
trix model? We believe that the answer is negative. Re- 
capitulating the sequence of arguments, it is evident that 
everything hinges on the absence of singularities outside 
the complex unit circle defined through the compactified 
energy argument. This in turn is a guaranteed feature as 
long as the single particle retarded SCBA Green function 
of the model system has a well defined pole structure be- 
low the real axis. In practice, for condensed matter sys- 
tems with non-universal high energy behavior, it may be 



difficult to find closed solutions of the mean field equa- 
tions that manifestly display this feature. We believe, 
however, that this is a practical rather than a principle 
difficulty. Summarizing, the main goal of part II was 
to show that the causality feature underlying this and 
previous analyses of the effective Keldysh action can be 
made an explicit ingredient of the model, on the expense 
of including large energy asympotics. In practical appli- 
cations of the formalism, one will normally use causality 
feature pragmatically, as exemplified in part I. Yet it can- 
not be excluded that situations arise, where large scale 
spectral structures become essential. 

To conclude, on the example of level statistics, we 
have demonstrated how non-perturbative quantum ef- 
fects may be incorporated into the framework of the dy- 
namic Keldysh cr-model. In many respects, the Keldish 
scheme appears to be simpler and more transparent than 
its relatives, replica and SUSY. We expect our methods 
to be useful in the analysis of interaction phenomena in 
disordered electronic systems. 

We benefited a lot from numerous discussions with A. 
Andreev and M. Janssen. A.K. was partially supported 
by the BSF-9800338 grant. 
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